Quantum transport through STM-lifted single PTCDA molecules 



00 

o 
o 

as: 
< 

OS 



Florian Pump, 1 Ruslan Tcmirov. 2 Olga Neucheva, 2 Scrguci Soubatch, 2 
Stefan Tautz, 2 Michael Rohlfing, 3 and Gianaurelio Cunibcrti 1 

'institute for Materials Science and Max Bergmann Center of Biomaterials, 
Dresden University of Technology, D-01062 Dresden, Germany 
2 Institut fur Bio-und Nanosysteme 3, JARA, Forschungszentrum Jiilich, D-,52425 Jiilich, Germany 
3 Department of Physics, University of Osnabriick, D-49069 Osnabriick, Germany 

(Dated: August 19, 2008) 

Using a scanning tunneling microscope we have measured the quantum conductance through 
a PTCDA molecule for different configurations of the tip-molecule-surface junction. A peculiar 
conductance resonance arises at the Fermi level for certain tip to surface distances. We have relaxed 
the molecular junction coordinates and calculated transport by means of the Landauer/Keldysh 
approach. The zero bias transmission calculated for fixed tip positions in lateral dimensions but 
different tip substrate distances show a clear shift and sharpening of the molecular chemisorption 
level on increasing the STM-surface distance, in agreement with experiment. 

PACS numbers: 68.37.Ef, 81.07.Pr, 73.63.-b, 68.43.Bc, 85.85.+J 
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I. INTRODUCTION 

Since the emergence of molecular electronics in the late 
nineties of the previous century enormous progress has 
been achieved in both experimental and theoretical un- 
derstanding of the underlying microscopic mechanisms 
governing the charge (and also recently spin) migration 
at the molecular scale PJ. In this sense, the realization 
of electronic functions at nanometer scales, as pioneered 
by the visionary work of Aviram and Ratner in 1973 Q, 
is slowly turning from a dream into a reality. However, 
there are still formidable challenges to be overcome. 

The current state of affairs can be characterized by 
the following observations: in nanometer-sized junctions 
the connection between the molecule and the electrode 
greatly affects the current- voltage characteristics [U , but 
there are currently no robust methods to image and de- 
termine the precise adsorption site and conformation of 
the molecule on this length scale Q. We are thus fac- 
ing a dilemma: on the one hand, we know from theory 
that contacts are extremely important H, @, 0] , on the 
other hand there is - within conventional approaches - 
no way to obtain atomic-scale information on these con- 
tacts. The statistical analysis of large data sets (thou- 
sands of experiments) has proven to be one way out of 
the dilemma. However, this procedure is evidently very 
laborious. For this reason, we wish to explore an alterna- 
tive approach in the present paper, namely to chemically 
contact single molecules in a reproducible way, such that 
reproducible transport spectra can be obtained even in 
small data sets. This approach is currently being ex- 
plored; see e.g. refs. [1, Q. It allows for transport ex- 
periments in structurally well-controlled junctions and 
thus offers an excellent interface to transport calcula- 
tions. Contacts to single atoms have also been stud- 
ied [II [II [11 03. 

Experimentally, there are three main techniques which 
form the pillars of the current state-of-the-art: mechan- 
ically controllable break junctions (MCBJ) 0, [H, [H, 




[l9l . I20I l2ll I22I [23j], electromigration techniques 
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Of course, each of these techniques has its specific ad- 
vantages and drawbacks. For instance, due to their rigid 
configuration, MCBJs are very stable from a mechanical 
point of view; consequently, experiments with variable 
temperatures can be performed. However, MCBJs have 
no lateral scanning ability and it is very difficult to get 
information on the detailed microscopic structure of the 
contacts (see above). Electromigration techniques allow 
for building three-terminal junctions in a straightforward 
way, but again the interface to the molecules is not very 
well defined. 

In contrast to the above techniques, STM-based trans- 
port experiments in the tip-molecule-surface geometry - 
by virtue of the STM imaging capability - allow for ex- 
perimental access to the atomic junction structure. How- 
ever, STM junctions arc mechanically less stable and, if 
operated in the tunnelling regime, they yield intrinsically 
asymmetric molecular junctions. Moreover, in an STM 




FIG. 1: Relaxed geometries for two different tip substrate 
distances cZts (A and B); top view of a PTCDA molecule 
(C). 
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junction electrical gating is often not available, because 
it is very difficult to realize a third electrode in the vicin- 
ity of tip and sample. Finally temperature dependent 
transport experiments are more difficult because of ther- 
mal drifts of the STM tip with respect to the sample 
surface. These problems notwithstanding, its is believed 
that STM-bascd transport experiments arc a worthwhile 
alternative to MCBJ and EM experiments, because they 
offer unique advantages, most importantly the structural 
control. As such, STM-based transport experiments are 
complementary to the other two techniques. Indeed, 
the STM tip cannot only address individual molecules, 
but also s pec ific sub-molecular groups or even individual 
atoms @, Hi EH, S3, S3, H, SH|. Moreover, it must be 
noted here that some of the disadvantages of the STM 
approach can in fact be overcome by using the STM in a 
slightly unconventional fashion: firstly, the contact asym- 
metry can be surmounted by contacting the molecule 
with the tip; this contact may either lead to the forma- 
tion of a strong covalent chemical bond which survives tip 
retraction ( "chemical contact" ) [|[ or be purely "mechan- 
ical" in the sense that on tip retraction the contact be- 
tween tip and molecule is opened again 0]. Secondly, in 
the case of chemical contacts the tip can also be used for 
mechanical (in contrast to capacitive) gating by chang- 
ing the effective coupling strengths between molecule and 
leads through variations of the vertical tip position Q. 
This might potentially allow one to scan different physi- 
cal regimes ranging from weak coupling (Coulomb block- 
ade), through intermediate coupling (Kondo physics) up 
to strong coupling (coherent transport). By its capabil- 
ity of manipulation, STM thus allows for tuning of the 
transport characteristics of individual molecules. 



As a way out of the dilemma concerning reproducible 
molecular electronics experiments mentioned above, we 
use here a two-step approach based on commensurate, 
highly ordered molecular layers on metal surfaces [1, H(| . 
In step 1, the powerful armoury of surface science is em- 
ployed to characterize the structural and electronic prop- 
erties of the metal-molecule contact. Once step 1 has 
been completed, the tip of a low-temperature STM is 
used to establish a covalent contact with an individual 
molecule of the highly ordered monolayer (step 2). Be- 
cause of the excellent imaging properties of the STM, it 
is in fact possible to select the part of the molecule which 
is contacted with very high accuracy. In this way a struc- 
turally well-defined molecular wire is realized which can 
be characterized with respect to its transport properties. 
In favorable cases, it may even be possible to apply a 
mechanical gating of the wire by tip movements @]. 

This paper is organized as follows: Sects. HT1 presents 
the experimental strategy and the most important results 
of a proof-of-principlc experiment; Sect. IIIIIIVI illustrate 
how the interface from these experiments to ab initio the- 
ory can be devised. 
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FIG. 2: Schematic representation of the bonding interaction 
of PTCDA with Ag(lll) [HEi]. 



II. STM EXPERIMENTS 

We have contacted a PTCDA (4,9,10-perylene- 
tetracarboxylic-dianhydridc) molecule with an STM tip, 
peeling it off the Ag(lll) surface and measuring its trans- 
port properties in the process. PTCDA/Ag(lll) is a 
prototypical model interface, for which the following is- 
sues have been addressed by previous work [47|: (1) 
The precise geometric and electronic structure of the ad- 
sorbed monolayer, including site-specific effects and influ- 
ence of the intermolecular interactions, have been deter- 
mined (cf. the bonding model in Fig. [5]). (2) The mech- 
anisms of the chemical bonding of PTCDA to Ag(lll), 
including orbital hybridization, charge transfer, local ver- 
sus extended bonds, and internal molecular distortion 
have been analyzed. (3) The image formation in the 
STM has been simulated from first principles. (4) Vari- 
ous density functionals, including GGA and LDA, have 
been evaluated for the PTCDA/Ag(lll) adsorption sys- 
tem. In particular, it must be noted here that according 
to the bonding model of Fig. [2] the carboxylic oxygen 
atoms are involved in local Ag-0 bonds. By their bond 
length they effectively define the vertical distance of the 
perylene core above the metal surface (50| . This in turn 
controls the degree of metal-molecule hybridization and 
charge transfer, because both proceed via the extended 
7r-electron system of the perylene backbone. At the same 
time, it has been shown that the Ag-0 bonds have only a 
secondary influence on the electronic properties of the ad- 
sorbed molecule @. To first order, we can thus conclude 
that the Ag-0 bonds act as (by themselves electronically 
inactive) mechanical clamps which control the electron 
interaction via the -as we will see below mechanically 
weaker - perylene- metal bond. 

PTCDA/Ag(lll) is a system suited well to test the 
two-step approach outlined at the end of the introduc- 
tion, because of various reasions. Firstly, the bond of 
PTCDA to the Ag(lll) substrate is very well character- 
ized, i.e. a lot of experiments with relevance to step 1 have 
already been performed. Secondly, because of the coex- 
istence of an electronically active but mechanically weak 
extended bond and an electronically inactive but mechani- 
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FIG. 3: Approach and retraction spectra of the conductance 
of the stretched tip-PTCDA-Ag(lll) (a, b) or tip-PTCDA- 
Au(lll)(c) transport junctions. Bias voltage Vb — 2 mV for 
(a, b) and 4 mV for (c). Images in (a-c) show the investigated 
molecule before and after the experiment, the cross indicates 
the point of approach. The zero of the z scale in (a-c) was 
established according to the procedure described in Ref. 0. 
Vertical bars in (a, b): range of differential conductance values 
from spectra in (d, e)). (d) Series of G = dl/dV vs. Vb) 
spectra (smoothed) measured at z positions indicated in (a) 
by vertical bars (V mo d = 4 mV, Vmod = 723 Hz), (e) Series 
of G — dl/dV vs. Vb spectra (smoothed) measured at z 
positions indicated in (b) by vertical bars (Vmod = 4 mV, 
fmod = 723 Hz). Colours in (d, e) code the regimes in (a, b). 

cally strong local bond (see Fig. [5]) , the molecule offers the 
possibility for docking the STM tip to the mechanically 
active terminal of the molecule and gradually cleaving the 
mechanically weaker, electronic bond to the substrate by 
tip retraction from the surface, thereby tuning the elec- 
tronic coupling strength of the molecule to the substrate. 
We shall see that in effect this indeed amounts to a me- 
chanical gating of transport through the molecular wire. 
In this simple way a single molecule transport experiment 
which combines a very high degree of structural control 
at the contacts with electronic tunability of one of the 
contacts is devised. 

First it must be established that it is indeed possible 
to contact the PTCDA molecule at one of the carboxylic 
oxygens. This can be concluded from the z-approach 
spectra, in which the conductance is measured as a func- 
tion of tip distance from the molecule Q. If the tip is 
approached above one of the carboxylic oxygen atoms, a 
sudden increase of conductance by more than one order of 
magnitude is observed at a tip distance of approximately 
2 A above the oxygen atom. This step in the conduc- 
tance is beyond the z-resolution of the STM. Since it is 
only observed above the carboxylic oxygen atoms, we can 
conclude that at the conductance step the oxygen atom 



flips up into contact with the tip, to a new potential 
minimum. Further analysis suggests that the tip-oxygen 
interaction can only be a covalent bond Q. For exam- 
ple, the reaction of the tip is independent of the bias 
direction, which excludes an electrostatic effect. Also, 
the tip-oxygen contact is strong enough to remove the 
molecules from the surface. Finally, after jumping into 
contact, the conductance reaches values which are typical 
of covalent bonds (in the range of 0.1 Go)[5l[ With these 
conductance values, currents through a single molecule 
can become as large as a few micro-amperes before the 
molecule is destroyed. 

Once the tip-oxygen contact has been formed, it is 
possible to retract the tip. We have mentioned already 
that in this way it is possible to remove the molecule 
from the surface completely. Analyzing the retraction 
more closely, we find that the conductance remains larger 
than the corresponding tunnelling values for retraction 
distances up to approximately 12 A Q. Hence we can 
conclude that in this range of tip positions the molecule 
is still present in the junction. Since the length of the 
molecule is 12 A, this indicates that the molecule is really 
lifted from a parallel to a vertical orientation by the tip, 
i.e. the junction can indeed be stretched out into a ver- 
tical wire. On tip retraction, two clearly distinguishable 
regimes are observed: During the first 4 A, the retraction 
spectrum is smooth, but in the range from 4 A to 12 A 
it becomes rather jerky. It is clear that if the molecule 
is lifted up into the vertical, its lower end must at some 
point slide over the surface, and this is not expected to 
be a smooth but rather a discontinuous movement. 

The possibility to achieve a vertical orientation of the 
molecule indicates that wc arc indeed able to modify the 
contact to the substrate in wide bounds. Most notably, 
for PTCDA on Ag(lll) we invariably observe that the 
conductance goes through a maximum within the first 
4 A of tip retraction. To understand the origin of this 
maximum, we have performed experiments in which tip 
retraction has been halted at regular intervals to measure 
full differential conductance spectra in the range from 
— 50 mV to 50 mV. An example of such a series of spec- 
tra is shown in Fig. [3^,, b, d, e. As the tip is retracted, 
we observe in the spectra the approach and sharpening 
of density of states from the left toward the Fermi level. 
(We note here that for not too strongly coupled wires, the 
energy-resolved transmission matrix and -accordingly - 
the differential conductance are proportional to the elec- 
tronic density of states in the wire). This density of 
states must be assigned to the former LUMO of free 
PTCDA, which serves as the bonding orbital responsi- 
ble for the electronic interaction (cf. the bonding sketch 
in Fig. [2]). Its behavior on junction stretching can be 
explained by reversed chemisorption, i.e. dehybridization 
of the LUMO from metal states (which leads to sharpen- 
ing) and an upward shift which follows the one-electron 
image potential of the occupied orbital (negative charge!) 
above the metal surface [|[ . Incidentally, this interpreta- 
tion is confirmed by analogous experiments with PTCDA 
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on Au(lll), which instead of the conductance maximum 
always show a monotonous decrease of junction conduc- 
tance on tip retraction, as shown in Fig. [3t. Since the 
interaction of PTCDA with Au(lll) is merely physisorp- 
tive [13, HH , little or no charge transfer between molecule 
and the Au(lll) substrate occurs, hence the LUMO of 
PTCDA is never pulled below the Fermi level. Lifting the 
molecule from the Au(lll) surface thus can only result 
in a further shift of the LUMO away from the Fermi level 
(due to decreased polarization screening). The absence 
of the maximum in the conductance curve for Au(lll) 
thus indicates that on Ag(lll) it must be linked to the 
electronic interaction of the molecule with the substrate, 
as described above. 

At the end of the tip retraction experiment, the density 
of states corresponding to the LUMO becomes pinned 
at the Fermi level Q. This cannot be explained in 
the framework of the reversed chemisorption model, for 
which one would naively expect the level to be pulled 
through the Fermi level, following the one-electron im- 
age potential further. At least it is difficult to ratio- 
nalize in a one-particle model why the level should be 
pinned exactly at the Fermi level, which corresponds to 
exactly one electron in the orbital. However, if many- 
body physics is taken into account, there is a mechanism 
that would naturally explain the observed pinning. This 
is the Kondo effect. In Ref. [|[ we have analyzed the pos- 
sibility to explain the pinning by the Kondo effect. To 
this end, an exact solution of the Anderson Hamiltonian 
541. using the numerical renormalization group (NRG) 
55, [HH, [13], has been used. Concerning the choice of 
the Anderson Hamiltonian as a model for the present 
system, we recall that in the formulation by Newns (48| 
this Hamiltonian is widely used as a model for molecular 
chemisorption on metal surfaces. Indeed it is known from 
our earlier work that for PTCDA/Ag(lll)-chemisorption 
the LUMO of PTCDA behaves essentially as a local level 
in the Newns- Anderson model 0, [13, HH . Accordingly, 
the Anderson Hamiltonian is an appropriate model for 
PTCDA/Ag(lll), if the role of many-particle correla- 
tions is to be analyzed. 

The NRG calculation shows that pinning of the bond- 
ing orbital at the Fermi level may indeed be observed 
even outside the fully developed Kondo regime, i.e. at the 
border to the mixed valence regime [H . This is an impor- 
tant finding, because our data suggest that on tip retrac- 
tion the tip/PTCDA/Ag(lll) wire remains on the bor- 
derline between the Kondo and mixed valence regimes. 
This can be concluded from the fact that we do not ob- 
serve a very strong sharpening of the zero bias feature 
with respect to the single particle peak. (The observed 
sharpening amounts to a factor of approximately 2 only.) 
As a consequence, we also do not resolve the two single 
particle peaks at eo (singly occupied LUMO) and eo + U 
(doubly occupied LUMO) separately, but rather a com- 
posite peak which consists of the moderately sharpened 
many body peak in the center and the two single par- 
ticle peaks in the left and right tails. In this context 



it is the most significant finding of the NRG solution 
that a pinning of the composite resonance to the Fermi 
level is indeed observed in the borderline regime between 
mixed valence and Kondo physics - in a parameter range 
that is suggested by our experiments. The experiment 
in Fig. [3] shows that before entering the pinning regime 
eo roughly amounts to 40mcV. Note that the order of 
magnitude of the intra-orbital Coulomb repulsion which 
we have employed (100 meV) has been corroborated by a 
many-body calculation using the GW formalism, which 
yields U = 0.25 eV for PTCDA on the silver substrate. 
This value includes a large contribution of the screening 
by the metal surface: for the free molecule a U of 3.5 eV 
has been calculated. 

In conclusion, the Kondo scenario in our experiments 
is supported by strong circumstantial evidence. It must 
be stressed, however, that an experimental test of great 
significance cannot be performed, since we cannot stabi- 
lize the stretched junction at elevated temperatures for 
times long enough to measure the conductance peak at 
zero bias. This means that the expected decay of the 
Kondo peak at temperatures around the Kondo tempera- 
ture (approximately 120 K in our case) cannot be verified 
Q. However, we note here that in an experiment with 
Ceo in a mechanically controlled break junction [5^ . in 
which the Au contacts transfer charge into the molecule 
and make the LUMO half-filled, a very similar zero bias 
conductance peak has been observed on opening the junc- 
tion, and in that case it could be verified that the peak 
indeed decays following the dependence predicted for a 
Kondo peak. Because the transport physics in the exper- 
iment by Parks et al. is very similar to the experiment on 
PTCDA/Ag(lll) reported here, it seems permissible to 
invoke an analogous temperature behavior for our junc- 
tions, too. Incidentally, in the experiments by Parks et 
al. only the pinning regime was observed. The cross-over 
from the single-particle to the many-particle regime was 
not reported there. 

In spite of the internal consistency of the Kondo sce- 
nario outlined above, a final proof is needed. Because 
of its structurally well-defined nature, the present sys- 
tem offers the unique possibility to provide this proof by 
a comparison between theory and experiment: Because 
experimentally the structure of the PTCDA/Ag(lll) 
adsorbate-substrate complex is well known, and because 
its structural and electronic properties have reliably been 
calculated by DFT, we can perform transport calcula- 
tions of the realistic junction. If this is done at vari- 
ous levels of sophistication, the physics of the pinning 
can be illuminated: if a single-particle calculation does 
not show pinning, it is very likely that the pinning is 
a many-body phenomenon. Conversely, if pinning is al- 
ready observed at the single-particle level, the assump- 
tion of Kondo physics may not be necessary, although the 
mechanism outlined above should still remain valid; its 
contribution to the observable behavior could then just 
be dominated by some other mechanism. 

In order to realize this program, two types of calcula- 
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tions are necessary: First, starting from the experimen- 
tally verified structure of the relaxed junction (no tip 
retraction), the structural and electronic evolution of the 
junction on tip retraction must be calculated. Secondly, 
on the basis of the equilibrium structures of the stretched 
junction, transport calculations may be performed. In 
the next section we shall turn to a DFT simulation of 
the process of lifting the molecule off the surface with the 
retracting STM tip. The results of this simulation will 
fully corroborate the qualitative picture with respect to 
contact formation and peel-off deduced from the experi- 
ments. 



III. JUNCTION ATOMISTIC CONFIGURATION 

The fundamental questions to be answered by the DFT 
simulation of the junction stretching process are the fol- 
lowing. Firstly: can we verify the flipping up of the car- 
boxylic oxygen atom toward the tip and the formation of 
a chemical bond between the two? Secondly: what hap- 
pens to the molecule structurally during tip retraction? 

The geometrical configuration of the present system, 
consisting of a metal substrate, molecule, and STM tip, 
can successfully be described using density-functional 
theory (DFT). DFT, as realized by the SIESTA package 
[60l . l6ll . !g3| , has already been used by some of the present 
authors to investigate the adsorption of a monolayer of 
PTCDA on the Ag(lll) surface 51, HI . It has turned 
out that PTCDA adsorbs on bridge sites above the sub- 
strate, at an equilibrium height of about 2.8 A, charac- 
terized by a nearly flat geometry with a slight molecular 
distortion. All these features are confirmed by the corre- 
sponding experimental data [H, [63| , including the par- 
tial occupation of the molecular LUMO state by electron 
charge transfer from the substrate. One should note that 
DFT as such does not appear to be fully appropriate for 
the adsorption of PTCDA on Ag(lll) which apart from 
the chemical interaction has a sizable contribution from 
physisorption, because long-range van der Waals corre- 
lation effects are absent from in DFT |49|. Nonetheless, 
the theoretical results obtained in Refs. [40l. l49l . [63| are 
in close agreement with the measured data, thus justify- 
ing the use of DFT for the questions to be studied here. 
This holds in particular for the local-density approxima- 
tion (LDA) to the exchange-correlation functional, which 
yields much better results than the generalized gradient 
approximation (GGA), which fails to bind the molecule 
to the surface. 

Based on the good results obtained previously for 
PTCDA on Ag(lll) [4!|, we here proceed to include an 
STM tip into the DFT calculation. Since we are not 
only interested in STM images (which can be described 
successfully in terms of a Tersoff-Hamann approach, i.e. 
without explicit inclusion of the tip [49l]). we focus on the 
mechanical interaction at close distances, i.e. we allow 
for the relaxation of the adsorbed molecules under the 
influence of forces exerted by the tip and vice versa. To 



this end, the tip is included as a pyramidal structure of 
ten tungsten atoms. In close relation to the experimental 
setup described in Ref. [|[ , we place the tip above the cor- 
ner of the molecule, with the tip apex atom exactly above 
one of the carboxylic oxygen atoms. A chemical bond be- 
tween the carboxylic oxygen atom and the apex tungsten 
atom forms, with a bond length of 2.1 A. The calcula- 
tions thus fully confirm the conclusion drawn previously 
from our experimental data. Thereafter, the STM tip 
is lifted up in steps of 0.2 A. After each movement of 
the tip, the molecule is allowed to relax and follow the 
tip. One should note that due to the large complexity of 
the molecule it requires up to 100 geometric relaxation 
steps until the molecule has adjusted to the retreated 
tip. Fig. [T] shows two typical relaxed structures of this 
procedure. In configuration A, the molecule is still close 
to its free adsorption structure. The carboxylic oxygen 
atom at the left-hand side of the molecule, however, has 
formed the chemical bond to the tungsten tip, which has 
a position of 5.4 A above the Ag surface, corresponding 
to the initial phase of tip retraction. Upon retreating 
the tip further from the surface, the carboxylic oxygen 
atom remains chemically bonded to the tip and follows 
it, causing the molecule to bend and peel off the surface. 
Fig. [U B displays the configuration after retreating the 
tip by 3.6 A. 

The simulation reveals that the junction structure is 
controlled by three competing energies and interactions, 
(i) The bond between the carboxylic oxygen and the tip 
apex atom appears to be quite stiff and to be of a covalent 
nature, (ii) The molecule is still bonded to the substrate 
by a delocalized interaction mechanism (cf. Fig. [5]). This 
mechanism is driven by the electron charge transfer from 
the substrate into the molecular LUMO state, which is 
delocalized over the entire perylene core of the molecule 
[49l ] . This causes attraction between every part of the 
molecule and the substrate, (iii) The perylene core is 
relatively soft and can easily be deformed. The sum of 
all three interaction mechanisms finally yields the curved 
structure of the molecule in Fig. [TJ with its right-hand 
part still arranged flat on the substrate and the left-hand 
part lifted off. Apparently, the attraction of the right- 
hand part of the molecule to the surface is still stronger 
than the energy cost of bending the molecule. Again this 
is in good agreement with the behavior deduced from 
experiments. At present, the calculations have not yet 
reached the regime in which the molecule slides across 
the surface. 

One important feature during the peeling process con- 
cerns the occupation of the LUMO state. This can be 
monitored by projecting the gas-phase LUMO wave func- 
tion onto the states of the substrate-adsorbate system. 
Results for three configurations (at tip-surface distances 
of 5.4 A, 6.8 A, and 9.0 A) are displayed in Fig. g] These 
are the configurations for which transport calculations 
will be discussed further below. It is clearly visible that, 
in accordance with the mechanism as described above, 
the LUMO is shifted to higher energy during the peeling 
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FIG. 4: Density of states ofthe LUMO state, projected onto 
the electronic states of the substrate-adsorbate complex, for 
configurations obtained at tip-surface distances of 5.4 A, 6.8 
A, and 9.0 A. 



process, and at 9.0 A tip height it has already a substan- 
tial density of states above the Fermi level. This is accom- 
panied by successive de-population and spectral sharp- 
ening. The completely removed molecule would have a 
sharp LUMO above the Fermi level of the substrate. 

The geometries as determined by the DFT simulation 
of the lift-off are the basis for the transport calculations 
described in the next section. 



IV. DFT-BASED COHERENT TRANSPORT 
CALCULATIONS 

As a further step towards a deeper understanding of 
the transport properties described in Sec. UH we perform 
transport simulations through a junction composed of 
substrate, molecule and STM tip. Even though model- 
based approaches often give valuable insight into under- 
lying transport mechanisms (64l . |6o1 [r36l . \6T\ , a realistic, 
material-specific description can only be achieved by an 
atomistic study. Given the size of the system, an ap- 
proximate density functional parametrized tight-binding 
method, DFTB, is used to account for the electronic 
structure of the hybrid system. The input geometries 
for these calculations are the fully relaxed ones obtained 
using the DFT package SIESTA as described above in 
Sec. IIIIl to obtain optimal reliability of the structure. 

As a first step in the transport simulations, we show in 
the present study as a characterizing feature the trans- 
mission function T (E) at zero bias voltage for differ- 
ent distances between the substrate and the STM tip. 
In this way one can follow the evolution of the molec- 
ular levels (peaks in the transmission spectra) and the 
electronic coupling to the contacts (metallic surface and 
STM tip, respectively) which is reflected by the width of 
these peaks. To do so, we use as computational method 
the quantum transport implementation gDFTB [68|, [6^] , 
which is a combination of the non-equilibrium Green 



function technique to describe the transport process 
and a DFT-based two-center charge self-consistent tight- 
binding method, DFTB [7(| , to account for the electronic 
structure of the extended molecular system (see Fig. [IJ . 
During the transport calculation, starting from the sys- 
tem Hamiltonian, a charge self-consistency loop over the 
non-equilibrium Green function, the density matrix and 
the real space charge density is performed under the non- 
equilibrium conditions imposed by the applied bias volt- 
age and the attached contacts. The resulting electric field 
is calculated via solving the Poisson equation at any loop 
cycle. The electric contacts are modeled as ideal semi- 
infinite wires, which are projected by decimation tech- 
niques onto an effective surface Green function. As the 
transport property considered here, the coherent trans- 
mission function is given by 



T {E ) = Tr [T 1 G T T 2 G a ] 



(1) 



where Ti and T 2 describe the coupling between the elec- 
trodes and the extended molecule at the contact layers. 
G T and G a are the retarded and advanced Green function 
of the extended molecule which are calculated as the re- 
sult of the self-consistency loop. 

It has to be stressed that the presented calculations are 
single-particle ones and cannot capture typical Kondo 
physics features or charging effects. However, the re- 
sults give deep insight into the electronic structure and 
its evolution during the lifting process and provide hints 
as to whether the observed experimental features can be 
explained in a single-particle picture or need a more so- 
phisticated explanation. 

In the transport calculations, a single PTCDA 
molecule, one surface layer consisting of 38 atoms and 
10 atoms of the top of the STM tip are taken as the 
extended molecule to be treated fully within DFTB for 
taking into account charge transfer processes between 
the surface and the molecule at all stages of the lifting 
process. Both leads are modeled by semi-infinite silver 
wires. It has to be noted that finite size effects may 
lead to a quantitative difference compared to studies of 
the system taking into account periodic boundary con- 
ditions which are currently carried out. The parameter 
sets used in this study were the ones available from the 
DFTB project [TO, [Hi, Ll2, Ll3]- Fig- El shows the zero 
bias transmission function T as a function of the incom- 
ing particle energy E for different tip-substrate distances. 
It can clearly be seen that the spectrum, especially the 
peak closest to the Fermi energy, moves in a reversed 
chemisorption process in a positive energy direction as 
observed in experiment. For the largest distance investi- 
gated in the present study, i.e. 9.0 A, the peak is located 
closely below the Fermi energy. The further evolution 
of the position of this peak with larger distances (the 
molecule lifted further from the surface) is of great inter- 
est as it could show whether the level crosses the Fermi 
energy or if it remains pinned there. The relaxation of 
the system for this parameter regime is currently under 
investigation. Furthermore it is evident from Fig. [5] that 
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cially of the frontier orbitals when lifting the molecule 
from the substrate. 

In continuation of the present work, it is necessary to 
explore stable system setups with larger tip-substrate dis- 
tances to follow the whole trace of the molecule during 
the lifting process as performed in the experiments, where 
also lateral movement of the molecules on the surface is 
observed. Furthermore, an investigation of the different 
behavior and properties comparing isolated molecules on 
a surface to a monolayer of PTCDA is presently per- 
formed. To finally clarify the origin of the experimentally 
observed features, it is necessary to find clear relations 
between the important magnitudes governing transport 
and the experimental input parameter, the tip substrate 
distance. In detail, these are the charging energy U, the 
coupling strengths between molecule and substrate, de- 
scribed by r, and the level energy e. Depending on their 
functional dependence during the lifting process, it will 
be possible to decide whether a single particle picture 
is sufficient to describe the observations or if a Kondo 
picture needs to be applied. 

Quite apart form the question whether the Kondo sce- 
nario is applicable for our model system, our proof-of- 
principle experiment and the simulation program, carried 
out here for PTCDA, provide a benchmark for the com- 
parison of a single-molecule transport experiment and 
ab-initio transport simulations, because hardly ever is 
structural information on the transport junction avail- 
able that is so detailed. 



FIG. 5: Transmission at zero bias as a function of energy for 
different tip substrate distances oIts- 



the width of the peaks in the spectra is decreased for 
larger distances, as, due to a reduction of the overlap of 
the molecular 7r-system with the metal electronic system, 
the electronic coupling strength and hence the injection 
rate between the molecule and the metallic substrate is 
reduced (Breit-Wigner physics). 



V. CONCLUSIONS AND OUTLOOK 

Using quantum transport methods, we have shown the 
motion of the energy levels in positive direction, espe- 
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